Computational approaches to investigate fluoride binding, selectivity and transport across the membrane.

The use of molecular dynamics (MD) simulations to study biomolecular systems has proven reliable in elucidating atomic-level details of structure and function. In this chapter, MD simulations were used to uncover new insights into two phylogenetically unrelated bacterial fluoride (F-) exporters: the CLCF F-/H+ antiporter and the Fluc F- channel. The CLCF antiporter, a member of the broader CLC family, has previously revealed unique stoichiometry, anion-coordinating residues, and the absence of an internal glutamate crucial for proton import in the CLCs. Through MD simulations enhanced with umbrella sampling, we provide insights into the energetics and mechanism of the CLCF transport process, including its selectivity for F- over HF. In contrast, the Fluc F- channel presents a novel architecture as a dual topology dimer, featuring two pores for F- export and a central non-transported sodium ion. Using computational electrophysiology, we simulate the electrochemical gradient necessary for F- export in Fluc and reveal details about the coordination and hydration of both F- and the central sodium ion. The procedures described here delineate the specifics of these advanced techniques and can also be adapted to investigate other membrane protein systems.

Fluoride Ion Binding and Translocation in the CLCF Fluoride/Proton Antiporter: Molecular Insights from Combined Quantum-Mechanical/Molecular-Mechanical Modeling.

CLCF fluoride/proton antiporters move fluoride ions out of bacterial cells, leading to fluoride resistance in these bacteria. However, many details about their operating mechanisms remain unclear. Here, we report a combined quantum-mechanical/molecular-mechanical (QM/MM) study of a CLCF homologue from Enterococci casseliflavus (Eca), in accord with the previously proposed windmill mechanism. Our multiscale modeling sheds light on two critical steps in the transport cycle: (i) the external gating residue E118 pushing a fluoride in the external binding site into the extracellular vestibule and (ii) an incoming fluoride reconquering the external binding site by forcing out E118. Both steps feature competitions for the external binding site between the negatively charged carboxylate of E118 and the fluoride. Remarkably, the displaced E118 by fluoride accepts a proton from the nearby R117, initiating the next transport cycle. We also demonstrate the importance of accurate quantum descriptions of fluoride solvation. Our results provide clues to the mysterious E318 residue near the central binding site, suggesting that the transport activities are unlikely to be disrupted by the glutamate interacting with a well-solvated fluoride at the central binding site. This differs significantly from the structurally similar CLC chloride/proton antiporters, where a fluoride trapped deep in the hydrophobic pore causes the transporter to be locked down. A free-energy barrier of 10-15 kcal/mol was estimated via umbrella sampling for a fluoride ion traveling through the pore to repopulate the external binding site.

Structural basis of pH-dependent activation in a CLC transporter.

CLCs are dimeric chloride channels and anion/proton exchangers that regulate processes such as muscle contraction and endo-lysosome acidification. Common gating controls their activity; its closure simultaneously silences both protomers, and its opening allows them to independently transport ions. Mutations affecting common gating in human CLCs cause dominant genetic disorders. The structural rearrangements underlying common gating are unknown. Here, using single-particle cryo-electron microscopy, we show that the prototypical Escherichia coli CLC-ec1 undergoes large-scale rearrangements in activating conditions. The slow, pH-dependent remodeling of the dimer interface leads to the concerted opening of the intracellular H+ pathways and is required for transport. The more frequent formation of short water wires in the open H+ pathway enables Cl- pore openings. Mutations at disease-causing sites favor CLC-ec1 activation and accelerate common gate opening in the human CLC-7 exchanger. We suggest that the pH activation mechanism of CLC-ec1 is related to the common gating of CLC-7.

Novel Mg 2+ binding sites in the cytoplasmic domain of the MgtE Mg 2+ channels revealed by X-ray crystal structures.

MgtE is a Mg 2+-selective channel regulated by the intracellular Mg 2+ concentration. MgtE family proteins are highly conserved in all domains of life and contribute to cellular Mg 2+ homeostasis. In humans, mutations in the SLC41 proteins, the eukaryotic counterparts of the bacterial MgtE, are known to be associated with various diseases. The first MgtE structure from a thermophilic bacterium, Thermus thermophilus, revealed that MgtE forms a homodimer consisting of transmembrane and cytoplasmic domains with a plug helix connecting the two and that the cytoplasmic domain possesses multiple Mg 2+ binding sites. Structural and electrophysiological analyses revealed that the dissociation of Mg 2+ ions from the cytoplasmic domain induces structural changes in the cytoplasmic domain, leading to channel opening. Thus, previous works showed the importance of MgtE cytoplasmic Mg 2+ binding sites. Nevertheless, due to the limited structural information on MgtE from different species, the conservation and diversity of the cytoplasmic Mg 2+ binding site in MgtE family proteins remain unclear. Here, we report crystal structures of the Mg 2+-bound MgtE cytoplasmic domains from two different bacterial species, Chryseobacterium hispalense and Clostridiales bacterium, and identify multiple Mg 2+ binding sites, including ones that were not observed in the previous MgtE structure. These structures reveal the conservation and diversity of the cytoplasmic Mg 2+ binding site in the MgtE family proteins.

ClC Chloride Channels in Gram-Negative Bacteria and Its Role in the Acid Resistance Systems.

Pathogenic bacteria that colonize the human intestinal tract have evolved strategies to overcome acidic conditions when they pass through the gastrointestinal tract. Amino acid-mediated acid resistance systems are effective survival strategies in a stomach that is full of amino acid substrate. The amino acid antiporter, amino acid decarboxylase, and ClC chloride antiporter are all engaged in these systems, and each one plays a role in protecting against or adapting to the acidic environment. The ClC chloride antiporter, a member of the ClC channel family, eliminates negatively charged intracellular chloride ions to avoid inner membrane hyperpolarization as an electrical shunt of the acid resistance system. In this review, we will discuss the structure and function of the prokaryotic ClC chloride antiporter of amino acid-mediated acid resistance system.

Uncovering the Mechanism of the Proton-Coupled Fluoride Transport in the CLCF Antiporter.

Fluoride is a natural antibiotic abundantly present in the environment and, in micromolar concentrations, is able to inhibit enzymes necessary for bacteria to survive. However, as is the case with many antibiotics, bacteria have evolved resistance methods, including through the use of recently discovered membrane proteins. One such protein is the CLCF F-/H+ antiporter protein, a member of the CLC superfamily of anion-transport proteins. Though previous studies have examined this F- transporter, many questions are still left unanswered. To reveal details of the transport mechanism used by CLCF, we have employed molecular dynamics simulations and umbrella sampling calculations. Our results have led to several discoveries, including the mechanism of proton import and how it is able to aid in the fluoride export. Additionally, we have determined the role of the previously identified residues Glu118, Glu318, Met79, and Tyr396. This work is among the first studies of the CLCF F-/H+ antiporter and is the first computational investigation to model the full transport process, proposing a mechanism which couples the F- export with the H+ import.

Small molecule inhibitors of intestinal epithelial anion exchanger SLC26A3 (DRA) with a luminal, extracellular site of action.

The anion exchanger protein SLC26A3 (down-regulated in adenoma, DRA) is expressed in the luminal membrane of intestinal epithelial cells in colon, where it facilitates the absorption of Cl- and oxalate. We previously identified a 4,8-dimethylcoumarin class of SLC26A3 inhibitors that act from the SLC26A3 cytoplasmic surface, and demonstrated their efficacy in mouse models of constipation and hyperoxaluria. Here, screening of 50,000 new compounds and 1740 chemical analogs of active compounds from the primary screen produced five novel classes of SLC26A3-selective inhibitors (1,3-dioxoisoindoline-amides; N-(5-sulfamoyl-1,3,4-thiadiazol-2-yl)acetamides; thiazolo-pyrimidin-5-ones; 3-carboxy-2-phenylbenzofurans and benzoxazin-4-ones) with IC50 down to 100 nM. Kinetic washout and onset of action studies revealed an extracellular site of action for the thiazolo-pyrimidin-5-one and 3-carboxy-2-phenylbenzofuran inhibitors. Molecular docking computations revealed putative binding sites for these inhibitors. In a loperamide model of constipation in mice, orally administered 7-(2-chloro-phenoxymethyl)-3-phenyl-thiazolo [3,2-a]pyrimidin-5-one (3a) significantly increased stool weight, pellet number and water content. SLC26A3 inhibitors with an extracellular site of action offer the possibility of creating non-absorbable, luminally acting inhibitors with minimal systemic exposure following oral administration. Our findings also suggest that inhibitors of related SLC26 anion transporters with an extracellular site of action might be identified for pharmacological modulation of selected epithelial ion transport processes.

Charge neutralization of the active site glutamates does not limit substrate binding and transport by small multidrug resistance transporter EmrE.

EmrE, a small multidrug resistance transporter from Escherichia coli, confers broad-spectrum resistance to polyaromatic cations and quaternary ammonium compounds. Previous transport assays demonstrate that EmrE transports a +1 and a +2 substrate with the same stoichiometry of two protons:one cationic substrate. This suggests that EmrE substrate binding capacity is limited to neutralization of the two essential glutamates, E14A and E14B (one from each subunit in the antiparallel homodimer), in the primary binding site. Here, we explicitly test this hypothesis, since EmrE has repeatedly broken expectations for membrane protein structure and transport mechanism. We previously showed that EmrE can bind a +1 cationic substrate and proton simultaneously, with cationic substrate strongly associated with one E14 residue, whereas the other remains accessible to bind and transport a proton. Here, we demonstrate that EmrE can bind a +2 cation substrate and a proton simultaneously using NMR pH titrations of EmrE saturated with divalent substrates, for a net +1 charge in the transport pore. Furthermore, we find that EmrE can alternate access and transport a +2 substrate and proton at the same time. Together, these results lead us to conclude that E14 charge neutralization does not limit the binding and transport capacity of EmrE.

Activating alternative transport modes in a multidrug resistance efflux pump to confer chemical susceptibility.

Small multidrug resistance (SMR) transporters contribute to antibiotic resistance through proton-coupled efflux of toxic compounds. Previous biophysical studies of the E. coli SMR transporter EmrE suggest that it should also be able to perform proton/toxin symport or uniport, leading to toxin susceptibility rather than resistance in vivo. Here we show EmrE does confer susceptibility to several previously uncharacterized small-molecule substrates in E. coli, including harmane. In vitro electrophysiology assays demonstrate that harmane binding triggers uncoupled proton flux through EmrE. Assays in E. coli are consistent with EmrE-mediated dissipation of the transmembrane pH gradient as the mechanism underlying the in vivo phenotype of harmane susceptibility. Furthermore, checkerboard assays show this alternative EmrE transport mode can synergize with some existing antibiotics, such as kanamycin. These results demonstrate that it is possible to not just inhibit multidrug efflux, but to activate alternative transport modes detrimental to bacteria, suggesting a strategy to address antibiotic resistance.

Integrated AlphaFold2 and DEER investigation of the conformational dynamics of a pH-dependent APC antiporter.

The Amino Acid-Polyamine-Organocation (APC) transporter GadC contributes to the survival of pathogenic bacteria under extreme acid stress by exchanging extracellular glutamate for intracellular γ-aminobutyric acid (GABA). Its structure, determined in an inward-facing conformation at alkaline pH, consists of the canonical LeuT-fold with a conserved five-helix inverted repeat, thereby resembling functionally divergent transporters such as the serotonin transporter SERT and the glucose-sodium symporter SGLT1. However, despite this structural similarity, it is unclear if the conformational dynamics of antiporters such as GadC follow the blueprint of these or other LeuT-fold transporters. Here, we used double electron-electron resonance (DEER) spectroscopy to monitor the conformational dynamics of GadC in lipid bilayers in response to acidification and substrate binding. To guide experimental design and facilitate the interpretation of the DEER data, we generated an ensemble of structural models in multiple conformations using a recently introduced modification of AlphaFold2 . Our experimental results reveal acid-induced conformational changes that dislodge the Cterminus from the permeation pathway coupled with rearrangement of helices that enables isomerization between inward- and outward-facing states. The substrate glutamate, but not GABA, modulates the dynamics of an extracellular thin gate without shifting the equilibrium between inward- and outward-facing conformations. In addition to introducing an integrated methodology for probing transporter conformational dynamics, the congruence of the DEER data with patterns of structural rearrangements deduced from ensembles of AlphaFold2 models illuminates the conformational cycle of GadC underpinning transport and exposes yet another example of the divergence between the dynamics of different families in the LeuT-fold.

Oligomeric state of the aspartate:alanine transporter from Tetragenococcus halophilus.

The aspartate:alanine exchanger family of membrane transporters includes industrially important transporters such as succinate exporter and glutamate exporter. No high-resolution structure is available from this family so far, and the transport mechanism of these transporters also remains unclear. In the present study, we focus on the oligomeric status of the aspartate:alanine antiporter (AspT) of Tetragenococcus halophilus, which is the prototype of this family. To investigate the oligomeric structure of AspT, we established a system that produces high yields of highly purified AspT and determined the oligomeric structure of AspT by analysis with size exclusion chromatography coupled with multi-angle light scattering and blue native PAGE and by comparison of the wild-type AspT with a single-cysteine mutant that forms spontaneous inter-molecular thiol crosslinking. All the results consistently support the notion that AspT is a homodimer in solutions and in membranes.

Principles of Alternating Access in LeuT-fold Transporters: Commonalities and Divergences.

Found in all domains of life, transporters belonging to the LeuT-fold class mediate the import and exchange of hydrophilic and charged compounds such as amino acids, metals, and sugar molecules. Nearly two decades of investigations on the eponymous bacterial transporter LeuT have yielded a library of high-resolution snapshots of its conformational cycle linked by solution-state experimental data obtained from multiple techniques. In parallel, its topology has been observed in symporters and antiporters characterized by a spectrum of substrate specificities and coupled to gradients of distinct ions. Here we review and compare mechanistic models of transport for LeuT, its well-studied homologs, as well as functionally distant members of the fold, emphasizing the commonalities and divergences in alternating access and the corresponding energy landscapes. Our integrated summary illustrates how fold conservation, a hallmark of the LeuT fold, coincides with divergent choreographies of alternating access that nevertheless capitalize on recurrent structural motifs. In addition, it highlights the knowledge gap that hinders the leveraging of the current body of research into detailed mechanisms of transport for this important class of membrane proteins.

Theoretical study of active secondary transport: Unexpected differences in molecular mechanisms for antiporters and symporters.

Successful functioning of biological cells relies on efficient translocation of different materials across cellular membranes. An important part of this transportation system is membrane channels that are known as antiporters and symporters. They exploit the energy stored as a trans-membrane gradient of one type of molecules to transport the other types of molecules against their gradients. For symporters, the directions of both fluxes for driving and driven species coincide, while for antiporters, the fluxes move in opposite directions. There are surprising experimental observations that despite differing only by the direction of transport fluxes, the molecular mechanisms of translocation adopted by antiporters and symporters seem to be drastically different. We present chemical-kinetic models to quantitatively investigate this phenomenon. Our theoretical approach allows us to explain why antiporters mostly utilize a single-site transportation when only one molecule of any type might be associated with the channel. At the same time, the transport in symporters requires two molecules of different types to be simultaneously associated with the channel. In addition, we investigate the kinetic constraints and efficiency of symporters and compare them with the same properties of antiporters. Our theoretical analysis clarifies some important physical-chemical features of cellular trans-membrane transport.

Structure, mechanism and lipid-mediated remodeling of the mammalian Na+/H+ exchanger NHA2.

The Na+/H+ exchanger SLC9B2, also known as NHA2, correlates with the long-sought-after Na+/Li+ exchanger linked to the pathogenesis of diabetes mellitus and essential hypertension in humans. Despite the functional importance of NHA2, structural information and the molecular basis for its ion-exchange mechanism have been lacking. Here we report the cryo-EM structures of bison NHA2 in detergent and in nanodiscs, at 3.0 and 3.5 Å resolution, respectively. The bison NHA2 structure, together with solid-state membrane-based electrophysiology, establishes the molecular basis for electroneutral ion exchange. NHA2 consists of 14 transmembrane (TM) segments, rather than the 13 TMs previously observed in mammalian Na+/H+ exchangers (NHEs) and related bacterial antiporters. The additional N-terminal helix in NHA2 forms a unique homodimer interface with a large intracellular gap between the protomers, which closes in the presence of phosphoinositol lipids. We propose that the additional N-terminal helix has evolved as a lipid-mediated remodeling switch for the regulation of NHA2 activity.

Isolation and functional characterization of a glucose-6-phosphate/phosphate translocator (IbG6PPT1) from sweet potato (Ipomoea batatas (L.) Lam.).

Sweet potato (Ipomoea batatas (L.) Lam.) is a good source of carbohydrates, an excellent raw material for starch-based industries, and a strong candidate for biofuel production due to its high starch content. However, the molecular basis of starch biosynthesis and accumulation in sweet potato is still insufficiently understood. Glucose-6-phosphate/phosphate translocators (GPTs) mediate the import of glucose-6-phosphate (Glc6P) into plastids for starch synthesis. Here, we report the isolation of a GPT-encoding gene, IbG6PPT1, from sweet potato and the identification of two additional IbG6PPT1 gene copies in the sweet potato genome. IbG6PPT1 encodes a chloroplast membrane-localized GPT belonging to the GPT1 group and highly expressed in storage root of sweet potato. Heterologous expression of IbG6PPT1 resulted in increased starch content in the leaves, root tips, and seeds and soluble sugar in seeds of Arabidopsis thaliana, but a reduction in soluble sugar in the leaves. These findings suggested that IbG6PPT1 might play a critical role in the distribution of carbon sources in source and sink and the accumulation of carbohydrates in storage tissues and would be a good candidate gene for controlling critical starch properties in sweet potato.

Molecular basis for redox control by the human cystine/glutamate antiporter system xc.

Cysteine plays an essential role in cellular redox homoeostasis as a key constituent of the tripeptide glutathione (GSH). A rate limiting step in cellular GSH synthesis is the availability of cysteine. However, circulating cysteine exists in the blood as the oxidised di-peptide cystine, requiring specialised transport systems for its import into the cell. System xc- is a dedicated cystine transporter, importing cystine in exchange for intracellular glutamate. To counteract elevated levels of reactive oxygen species in cancerous cells system xc- is frequently upregulated, making it an attractive target for anticancer therapies. However, the molecular basis for ligand recognition remains elusive, hampering efforts to specifically target this transport system. Here we present the cryo-EM structure of system xc- in both the apo and glutamate bound states. Structural comparisons reveal an allosteric mechanism for ligand discrimination, supported by molecular dynamics and cell-based assays, establishing a mechanism for cystine transport in human cells.

Simulation of pH-Dependent Conformational Transitions in Membrane Proteins: The CLC-ec1 Cl-/H+ Antiporter.

Intracellular transport of chloride by members of the CLC transporter family involves a coupled exchange between a Cl- anion and a proton (H+), which makes the transport function dependent on ambient pH. Transport activity peaks at pH 4.5 and stalls at neutral pH. However, a structure of the WT protein at acidic pH is not available, making it difficult to assess the global conformational rearrangements that support a pH-dependent gating mechanism. To enable modeling of the CLC-ec1 dimer at acidic pH, we have applied molecular dynamics simulations (MD) featuring a new force field modification scheme-termed an Equilibrium constant pH approach (ECpH). The ECpH method utilizes linear interpolation between the force field parameters of protonated and deprotonated states of titratable residues to achieve a representation of pH-dependence in a narrow range of physiological pH values. Simulations of the CLC-ec1 dimer at neutral and acidic pH comparing ECpH-MD to canonical MD, in which the pH-dependent protonation is represented by a binary scheme, substantiates the better agreement of the conformational changes and the final model with experimental data from NMR, cross-link and AFM studies, and reveals structural elements that support the gate-opening at pH 4.5, including the key glutamates Gluin and Gluex.

Molecular mechanism of proton-coupled ligand translocation by the bacterial efflux pump EmrE.

The current surge in bacterial multi-drug resistance (MDR) is one of the largest challenges to public health, threatening to render ineffective many therapies we rely on for treatment of serious infections. Understanding different factors that contribute to MDR is hence crucial from the global "one health" perspective. In this contribution, we focus on the prototypical broad-selectivity proton-coupled antiporter EmrE, one of the smallest known ligand transporters that confers resistance to aromatic cations in a number of clinically relevant species. As an asymmetric homodimer undergoing an "alternating access" protomer-swap conformational change, it serves as a model for the mechanistic understanding of more complex drug transporters. Here, we present a free energy and solvent accessibility analysis that indicates the presence of two complementary ligand translocation pathways that remain operative in a broad range of conditions. Our simulations show a previously undescribed desolvated apo state and anticorrelated accessibility in the ligand-bound state, explaining on a structural level why EmrE does not disrupt the pH gradient through futile proton transfer. By comparing the behavior of a number of model charged and/or aromatic ligands, we also explain the origin of selectivity of EmrE towards a broad class of aromatic cations. Finally, we explore unbiased pathways of ligand entry and exit to identify correlated structural changes implicated in ligand binding and release, as well as characterize key intermediates of occupancy changes.

Slipknotted and unknotted monovalent cation-proton antiporters evolved from a common ancestor.

While the slipknot topology in proteins has been known for over a decade, its evolutionary origin is still a mystery. We have identified a previously overlooked slipknot motif in a family of two-domain membrane transporters. Moreover, we found that these proteins are homologous to several families of unknotted membrane proteins. This allows us to directly investigate the evolution of the slipknot motif. Based on our comprehensive analysis of 17 distantly related protein families, we have found that slipknotted and unknotted proteins share a common structural motif. Furthermore, this motif is conserved on the sequential level as well. Our results suggest that, regardless of topology, the proteins we studied evolved from a common unknotted ancestor single domain protein. Our phylogenetic analysis suggests the presence of at least seven parallel evolutionary scenarios that led to the current diversity of proteins in question. The tools we have developed in the process can now be used to investigate the evolution of other repeated-domain proteins.

A solid-supported membrane electrophysiology assay for efficient characterization of ion-coupled transport.

Transport stoichiometry determination can provide great insight into the mechanism and function of ion-coupled transporters. Traditional reversal potential assays are a reliable, general method for determining the transport stoichiometry of ion-coupled transporters, but the time and material costs of this technique hinder investigations of transporter behavior under multiple experimental conditions. Solid-supported membrane electrophysiology (SSME) allows multiple recordings of liposomal or membrane samples adsorbed onto a sensor and is sensitive enough to detect transport currents from moderate-flux transporters that are inaccessible to traditional electrophysiology techniques. Here, we use SSME to develop a new method for measuring transport stoichiometry with greatly improved throughput. Using this technique, we were able to verify the recent report of a fixed 2:1 stoichiometry for the proton:guanidinium antiporter Gdx, reproduce the 1H+:2Cl- antiport stoichiometry of CLC-ec1, and confirm loose proton:nitrate coupling for CLC-ec1. Furthermore, we were able to demonstrate quantitative exchange of internal contents of liposomes adsorbed onto SSME sensors to allow multiple experimental conditions to be tested on a single sample. Our SSME method provides a fast, easy, general method for measuring transport stoichiometry, which will facilitate future mechanistic and functional studies of ion-coupled transporters.